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1. Introduction 


The ABS conducts a number of repeated business surveys. One of the key objectives of these 
business surveys is to produce reliable estimates of change from one time period to the next. 
Output editing is usually the final check for major errors in the survey data to ensure that 
the estimates are reliable. There are a number of tools that can be used to check whether the 
estimates of change appear correct, such as comparing the estimate of change with external 
sources, examining the distributions of the survey data, or decomposing the movement 
estimates into components which lead to a better understanding of the estimates of change. 


The estimates of change are driven by a number of different factors, including: 


¢ birthing of new units, deathing of old units and changing characteristics of 
continuing units in the population of interest; 


¢ rotation of new units into sample and rotation of previously selected units out of 
sample; 


¢ changes in the values of auxiliary benchmark variables at the unit and aggregate 
levels; and 


¢ changes in the values of the variable of interest at the unit level. 


Some of these factors arise from actual changes in the population of interest, while others 
arise from a changes in the relationships between the population and the chosen samples. 


The ABS regularly updates its business survey frames by adding births, removing deaths 
and updating stratification variables. These frame changes have a direct impact on the 
estimates of change, as well as an indirect impact through the rotation of units into and out 
of the sample. The rotation of units into and out of the sample is also be caused by the ABS 
policy to rotate all small and medium business out of individual surveys after three years. 
This policy is upheld using synchronised sampling, a form Permanent Random Number 
(PRN) sampling procedure, developed by the ABS to minimise or maximise overlap 
between surveys as well as control rotation within surveys (Brewer, Gross and Lee 1999). 
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The ABS has recently developed a Generalised Estimation System for processing its large 
scale annual and sub-annual business surveys, which allows the specification of an 
estimator from a wide group of estimators produced under generalised linear regression 
models. For those business surveys which produce estimates using auxiliary information in 
the form of known auxiliary totals, changes in the values of auxiliary benchmark variables at 
the unit and aggregate levels will have an impact on the estimates of change. 


In order to better understand the influence of these factors on the estimates of change, a 
methodology has been developed to decompose the movement estimates into a number of 
factors. Holt and Skinner (1989) present a framework for the decomposition of the net 
difference in the population means into effects due to net differences of means within 
domains and effects due to domain composition. The framework is also extended to allow 
for non-stationary populations with birthing of new units, deathing of old units and 
changing characteristics of continuing units. This paper presents an more comprehensive 
alternative methodology for the decomposition of the movement estimates which can be 
used as a diagnostic tool to assist in the output editing of sample surveys. 


Section 2 introduces the generalised regression estimator used in ABS business surveys. 
Section 3 presents one possible decomposition of the movement estimates into primary and 
secondary factors. Section 4 describes how the bootstrap variance estimator can be used to 
produce approximate variances for the various decomposition components. Section 5 
measures the biases and mean squared errors of the decomposition components and the 
bootstrap variance estimator of the decomposition components in a simulation study. 
Section 6 examines the ability of the proposed decomposition methodology to identify 
irregularities. 


QUESTION 1: Is the proposed decomposition methodology into primary and secondary 
factors the “best” approach to the decomposition of the movement estimates? 


ps Estimation for ABS Business Surveys 
2.1 The Generalised Regression Estimator 


Consider a finite population U divided into H strata U = {U1,..., Uj, ..., Un}, where U,;, is 
comprised of N;, units. The objective is to estimate the population total Y =. %, Yn, where 
yni is the value of the variable of interest y for unit i in stratum h. Assume there exists a set of 
auxiliary variables x, = {X1i,..., Xki, .., Xxi}] for which the population totals X= x are known. 
Suppose stratified simple random samples {s1,..., Si, ..., Su} without replacement (SRSWOR) 
of sizes {11,...,n,...,H} are drawn with selection probabilities 2; = n,/N;, from populations 
{U1,..., Uk, ..., Uy}. The sampling weights w; = 1/7; are those used in the Horvitz-Thompson 
estimator Y =h > wiy1i for variable of interest y. The generalised regression estimator is 


given by (Sarndal, Swensson and Wretman, 1992): 


Y=E wyit(X-Zwix, )P=Z wey, (2.1) 
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where g; is the g-weight for unit i, defined as: 


and c; are specified positive factors that relate to the variance structure of the linear 
regression model associated with the GREG estimator: 


yiax, Bei 
where E(e;) = 0, Var(e;) = cia” and Cov(ei, ¢;) = 0, for all i #j. 
The generalised regression estimator for domain of interest c is given by: 


Y= d wigiyi (2.2) 


1€S¢ 
where s, is the sample of units within domain of interest c. 
2.2 The Generalised Regression Estimator for Movement Estimates 
Define the population at time t by U = {U!),..., U)?, ..., Uf?}, where U” is comprised of N}? 


units. The objective is to estimate the difference ny the popalnon totals between the two 
time points, Y") = Y? —y =z » ie = 2 Yn , where ye ) is the value of the variable of 


ic?) 
interest y at time point t for unit i in ae h. The generalised regression movement 
estimator is given by: 
fm) @) (2), (2) (1) (0), (0) 
hg =, Wi SFY; a, W; Si; (2.3) 
The generalised regression movement estimator for domain of interest c is given by: 


on) (2)_(2), (2) (1) (0), (1) 
Ye =% wi'8i Yi —% wi Si Vi (2.4) 
1ES¢ 1€S¢ 


where s( is the sample of units at time point t within domain of interest c. 
2.3 The Movement Estimator as a Function of Four Factors 


The movement estimator (2.4) can be written as the difference between a product of four 
factors: 
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om) (2) (2) (2)—(2) DeaVsoo() 
Yo =UNG Pre She Vic 4 Ni Pre Sie Yc 2) 


where p is the estimated proportion of units at time point t in stratum h in domain of 
interest c: 


w 
a(t) _ ies\e 


he y w'?) 


a) is the weighted average g-weight at time point t in stratum h in domain of interest c: 


(t) , (t), (t) 
patti 
She = OO 

GP Pott 


and 9 is the average value of the variable of interest y at time point t in stratum h in 


domain of interest c: 


Y w? (t) 
_() ies - 


he ~ > w'?) 
1€S ie 

In this situation the average g-weight is weighted by the value of the variable of interest, 

which will mean that changes in the value of the variable of interest will impact on the 

average value of the variable of interest as well as the weighted average g-weight. The main 

reason for using the weighted average g-weight, rather than an unweighted average 

g-weight, is to ensure that the product of the four factors sum to the survey estimates. 

Alternatively, a. could be calculated as the average g-weight at time point t in stratum h in 

domain of interest c: 

oy w? ba 


ia (t) 
—(t) sexs, IES jyc 


he x w 


ies ) 


and i, could be calculated as the weighted average value of the variable of interest y at 
time point t in stratum h in domain of interest c: 
yr wi gy 


«__(t) 
—(t) _ ieshe 


Yn yp OOH 
x, a Si 


1€8 j, 
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which would also ensure that the product of the four factors sum to the survey estimates. In 
this situation the weighted average value of the variable of interest is weighted by the 
g-weight, which will mean that changes in the g-weight will impact on the average g-weight 
as well as the weighted average value of the variable of interest. The weighted average 
g-weight is used in the rest of the paper. It was chosen over the weighted average value of 
the variable of interest because changes in the estimates are more often driven by changes in 
the average value of the variable of interest. The results of the decomposition of the 
movement estimates will be easier to interpret if they are based on the unweighted average 
value of the variable of interest. 


QUESTION 2: Is appropriate to weight the average g-weight by the value of the variable 
of interest or weight the average value of the variable of the interest by the g-weight to 
ensure the product of the four factors sum to the survey estimates? 


3. Decomposition of Movement Estimates for ABS Repeated Surveys 
3.1 Decomposition into Four Primary Factors 


Using a standard decomposition methodology the movement estimator (2.5) can be 
decomposed into the following terms: 


a (m) a(1) -(2) (1)—(1) 
us =2 [NY 4 NY 1p Ps, cm i. + aa py ’|N Ni Bic Vic 


=( a (1)—(1) -(1)—(1) 
2 [zn i IN NIP Tn. +h in - -Fj, IN NP Be She ; (3.1) 


a(2) al) )—(1) =(2 a(1)—(1) 
+h [Ni —Ni? 1p, —Pi. Iss yi, +h [Ni -NP][3? - -3, aie U, 


—( =(1) Ja(i—( a ( ~(1) (2) =( 1)=(1) 
+h INP -NP Ione - Tn. Pn. By. +E x [pre Pre Ess — By IN Ni Uh 


4E [Pre ~Bhe Fie ~Tae INV Bic +E [Bie ~ Bic ke ~The PicVne -2) 


+EINP -Ni1 [927-3 
FE INY - NIP UP ie — Pre WGhe ~ Hie Bre 
+h [Ni = Ny ze -3n li -T, Pe 
45 Be —Bhe Bie ~Bhe Ge ~Tre INV (3.3) 
FEINY -NP UP ie ~Phe IBhe ~Zie Vie —Fre | (3.4) 
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The first four terms (3.1) are the main effects for the four factors, while the other terms (3.2), 
(3.3) and (3.4) are the two, three and four way interaction terms between the four factors. 
These interaction terms are difficult to interpret and are usually relative small compared to 
the main effects (Holt and Skinner, 1989). Although this decomposition methodology has 
the primary factors standardised at the first time period, it would have also been acceptable 
to have the primary factors standardised at the second time period. 


In order to avoid including the interaction terms, Das Gupta (1991) developed an alternative 
decomposition methodology, along the lines suggested by Kitagawa (1955). Using the 
alternative decomposition methodology the movement estimator (2.5) can be decomposed 
into the following four primary effects: 


Ye =E [NY -NY19(Bic Phe’ Bie Bic Vie Vic) (3.5) 
+3 [De Be JoNP NP Be Bic Tne Tre) (3.6) 
+E [Fi Fie OND NP Bie Pac Tier Tne) (3.7) 
+h li? - -7. Jo Ni? ,N NO iPS 8.) (3.8) 


where gla”, a), po b2, co) c)) = BaP bY + qu pV ee 4 qUp@eM 4 qWp2) 2 


a bY) 4 g@pUe® 4 q@@p@%eW 4 3q@p@e@y12 


The first terms (3.5) measures the effect of changes in the number of units on the survey 
frame, the second term (3.6) measures the effect of changes in the estimated proportion of 
units in the domain of interest, the third term (3.7) measures the effect of changes in the 
weighted average g-weight, and the fourth term (3.8) measures the effect of changes in the 
average value of the variable of interest. 


3.1.1. Empirical Evaluation of the Contribution of Interaction Terms 


The Retail Business Survey was used to evaluate the relative size of the interaction terms to 
determine whether is it appropriate to distribute these interaction terms into the main 
effects. The Retail Business Survey produces monthly estimates of the value of turnover for 
retail businesses classified by state and industry. The principal objective of the Retail 
Business Survey is to show month to month movements of turnover for retail industries. 


The contribution of the main effects and interaction terms to the relative movements 
between the June and July 2005 estimates of retail turnover are presented in Table 1, while a 
comparison of the main effects and the primary effects to the relative movements between 
the June and July 2005 estimates of retail turnover are presented in Table 2 
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Table 1: Contribution of Main Effects and Interaction Terms to Relative Movements‘ 


Net Contribution Percentage Contribution 
to Relative Movements to Relative Movements 
State 2nd Order | 3rd Order 2nd Order | 3rd Order 
Main Effects} Interaction | Interaction Relative |Main Effects| Interaction | Interaction 
Terms Terms Movements Terms Terms 
NSW 0.89 0.13 -0.05 0.97 91.4% 13.4% -4.8% 
Vic 0.54 0.25 -0.08 0.71 76.7% 34.6% -11.3% 
Qld 4.85 -0.11 0.00 4.73 102.4% -2.4% -0.0% 
SA 3.78 -0.25 0.02 3.55 106.6% -7.2% 0.5% 
WA 1.91 0.03 0.00 1.93 98.6% 1.4% 0.0% 
Tas 3.36 -0.05 0.19 3.51 96.0% -1.5% 5.5% 
NT 6.07 -0.46 0.00 5.62 108.1% -8.1% 0.0% 
ACT -0.16 0.93 -0.22 0.55 -29.0% 169.1% -40.1% 
Aust 1.99 0.08 -0.03 2.03 97.9% 3.8% -1.7% 


* All main effect and interaction terms which involve differences in the estimated proportion of units within the 


domain of interest (i.e. State) between the two time periods will be equal to zero, since the Retail Business Survey 
is stratified by State. 


Table 2: Comparison of Main Effects and Primary Effects 


Net Contribution Net Contribution 
to Relative Movements to Relative Movements 
of Main Effects of Primary Effects 
State Weight Average Weight Average 
Frame Change Value Frame Change Value 
Change Adjustment Change Change Adjustment Change 
Effect Effect Effect Effect Effect Effect 
NSW 0.81 0.38 -0.30 0.81 0.44 -0.27 
Vic 0.44 -0.30 0.41 0.48 -0.23 0.46 
Qld 0.34 2.22 2.28 0.31 2.20 2.23 
SA 0.72 2.62 0.45 0.58 2.58 0.38 
WA 0.79 1.11 0.01 0.76 1.14 0.03 
Tas 0.35 0.33 2.68 0.37 0.41 2.72 
NT 0.93 0.30 4.85 0.71 0.29 4.62 
ACT -0.01 3.18 -3.33 0.13 3.47 3.05 
Aust 0.59 0.87 0.53 0.58 0.91 0.54 


In most cases the percentage contribution to the relative movements of the main effects was 


between 90 and 110 percent. In those cases where the interaction terms contributed more 


than 10 percent of the relative movements, the relative movements were small (i.e. less than 


one percent) as were the net contribution to the relative movements. Furthermore, the net 


contribution to the primary effects are similar to the net contribution to the main effects. 


QUESTION 3: Is it appropriate to distribute the interaction terms into the main effects? 
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3,2 Decomposition into Secondary Factors 


The movement estimates are influenced by more factors than just the four primary factors 
described in Section 3.1. For example, stratum level changes in the number of units on the 
survey frame can occur as a births and deaths in the population as well as units changing 
stratum. Furthermore, changes in the average value of the variable of interest can occur as a 
result of changes in the values of the common sample units as well as differences between 
the values of the uncommon units. 


A pictorial illustration of the stratum level changes in the survey frame and survey sample 
between the two time periods is presented in Figure 1. 


Figure 1: Pictorial Illustration of Stratum Level Changes in the 
Survey Frame and Survey Sample 


Frame Time 1 


Frame 


Common Frame 
Death 


Sample Rotation Out 


Common Sample 


Frame Sample Rotation In 


Flip Out 


Frame Time 2 


The units on the survey frames within the stratum in the two time periods are represented 
by the left and right large rectangles respectively. The units on the common frame within the 
stratum are represented by the intersection of the two large rectangles, while the units not 
on the common frame within the stratum can be split into frame births, frame deaths, frame 
flips in and frame flips out. Those units which move into and out of the scope of a survey 
will be considered frame births and frame deaths. 


The units in the survey samples within the stratum in the two time periods are represented 
by the left and right ovals respectively. The units in the common survey sample within the 
stratum are represented by the intersection of the two ovals, while the units not in the 
common survey samples within the stratum can be split into sample rotations in, sample 
rotations out, sample flips in, sample flips out, sample births and sample deaths. 
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The purpose of the decomposition into secondary factors is to split the effect of changes in 
the number of units on the survey frame into frame births, frame deaths, frame flips in and 
frame flips out effects, and the effect of changes in the average value of the variable of 
interest into common sample, sample rotations in, sample rotations out, sample flips in, 
sample flips out, sample births and sample deaths effects. 


3.2.1 Frame Change Decomposition 


The difference in the number of units on the survey frames between the two time periods, 
can be written as: 


¢ differences in the number of frame births and frame deaths (i.e. frame growth); 
plus 


¢ differences in the number of frame flips in and frame flips out (i.e. frame stratum 
flips). 
(IN? ae Nw a [Ne _ NS 4s (ING Ne] 


Therefore, the frame change effect( 3.5) can be decomposed into a frame growth effect (3.9) 
and a frame stratum flip effect (3.10): 


(mf) (2) ()7,20) «(@) ~G) ~@) —G) —@) 
Y, =i [Nj —N;,, 1O(Pre+ Pre Sie She Ync Vine ) 


) (1d) a(1) a(2) =) =~) =) =(2) 
=> [NP —Nj"4 1P(P ic Phe She Bnet Vine Vie ) (3.9) 


i 0 a(1) -(2) <G) ~@2) —G) —() 
+h [Ni —Ni19( Bice Bc Sic She? Whe Vac ) (3.10) 


3.2.2 Average Value Change Decomposition 


The difference in the average value of units in the survey samples between the two time 
periods can be written as a function of: 


¢ differences in the average value of common sample units between the two time 
periods; 


¢ differences in the average value of sample rotations in and common sample units 
at the second time period; 


¢ differences in the average value of sample flips in and common sample units at 
the second time period; 


¢ differences in the average value of sample births and common sample units at the 
second time period; 


¢ differences in the average value of sample rotations out and common sample 
units at the first time period; 


¢ differences in the average value of sample flips out and common sample units at 
the first time period; and 
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¢ differences in the average value of sample deaths and common sample units at 
the first time period. 


[y= 7 |= rie q. | 


Vie he 
(2ri) (2fi) (2b) 
Nhe [—(ri) —(2c) Nhe (2fi) — —(2c) Nne [—(2b)  —(2c) 
as i he 4 he |+ 97 Vic |+ onl Vue | 
he he he 
(1ro) (1fo) ni d) 
Nhe [—(c) =(ro)] , Mne fac) (po) he (Ic) _ = (1d) 
+ —| he ~Yne ]+ (lo) li he Yn ]+- scone he | 
he Nic Nhe 


Therefore, the average value change effect (3.8) can be decomposed into a common sample 
unit effect (3.11), a sample rotations in effect (3.12), a sample stratum flip in effect (3.13), a 
sample birth effect (3.14), a sample rotations out effect (3.15), a sample stratum flip out effect 
(3.16) and a sample death effect (3.17): 


(mu) =) —(1) (1) ) aQ) aQ) —() —(2) 
Ye => | Yn lone, N Ny S Dict Pic/ rai ic) 
_y foc) =e) (1) NV (2) 2) a@) =) x~(2) 
=2 Tin Vac Jo(NPN h Prot Phe She? She ) (3.11) 
ng” ori) 00) (1) x72) al) a2) —G) =@) 
C —(2ri —(2c 1 2) 2a Z = a 
+h 2) li. —Vne Jo(NW?,N Nis Pres Phe Bie She ) (3.12) 
he 
z nj l7 (ff) = 10 NOON y (0) @)-= 4) =() 3.13 
as 12 Vic -Hin (Nin Nt Ba ay ora (3.13) 
Nic 
oe (2b) (2c) (1) x72) -@) ~2@) =) =@) 
C = —(2c a “a = = 
+h 12 li, Vote |p (Ni Ni Phe Pict Shc 7 he ) (3.14) 
he 
ae (1c) (1r0) ) a) ~Q@) =) =) 
C —(Ic ro N@ a a = 
+h ni La — Vue lo (Nj), Ni Pret Phe She She ) (3.15) 
he 
ny y (1c)  —(1fo) (1) x72) ~@) ~@) ~G) ~(@) 
+h nite) Fie Vne Jon} Nj Prot Pre She Sic) (3.16) 
ne (1c) (1d) (1) N?) (1) AQ) ~@) ~(2) 
C —Uc = 2) a as = = 
x ) ce Vn JocNi?, N h 1PhetPhe8he 1 &he ) (3.17) 


A description of the decomposition components (3.5 - 3.17) is presented in Appendix 1. 


Under this average value change decomposition, the common sample unit effect (3.11) could 
be considered to be inconsistent with the other average value effects (3.12 - 3.17). The 
common sample unit effect compares the differences in the average value of common 
sample units between the two time periods, while the other average value effects compare 
differences in the average value of uncommon sample units and common sample units at a 
single time period. 
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There are two main reason for adopting this average value change decomposition: (i) to 
ensure that the average value effects sum to the average value change effect (3.8); and (ii) it 
would appears to be logical to compare the average value of uncommon sample units with 
the average value of common sample units at a the same time period, because if the 
uncommon units have the same values as the common units then the decomposition effects 
would be equal to zero. 


In order for the other average value effects to be consistent with the common sample unit 
effect, then the other average value effects would need to compare differences in the average 
value of uncommon sample units at one time period and common sample units at the other 
time period. However, the average value effects would not sum to the average value change 
effect. 


QUESTION 4: Is it necessary for the common sample unit effect to be consistent with the 
other average value effects? 


In order to better interpret changes in the movement estimates, the primary and secondary 
factors can be collapsed into the following six broad components: 


(muc) xv(mbd) A(msf) Amr) amd) A(mw) 
+Y, +Y, +Y, +Y¥,. =+Y, (3.18) 


a (muc) a(mbd) (mfg) (ub) (mud) 
where y is the common sample unit effect (3.11), Y, ae Y, ‘s +Y, : +Y, is the net 
a(ms almfg) A (mufi) ~~ (ufo) 
birth-death effect (3.9, 3.14, 3.17), Y, Ve Y, e ne ae y +Y, i is the net stratum flip effect 
A(mr) A(muri) A (muro) a (md) 
(3.10, 3.13,3.16), Y, =Y, +Y;, is the net rotations effect (3.12,3.15), Y. is the 


a (mw) 


domain change effect (3.6), and Y, is the weight adjustment change effect (3.7). 
3.3 Decomposition of Movement Estimates of Rates 


The movement estimates of a rate between the two time periods can be written as a function 
of differences between the estimates of the numerator between the two time periods and 
differences between the estimates of the denominator between the two time periods. 
Therefore, the movement estimates of a rate can be decomposed into a numerator effect and 


a denominator effect: 


pilD) ~ (1) 
a(m) A. A. 
Yo = Bg” = Ea 
a(1l) a2) 
42) ~(1) (B. +B.’) 
= (a. al) — (3.19) 
Cc B, 
a(1) A(2) 
a2) =(1) (a. +A. 
+|-(B -B. er urcm (3.20) 
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The numerator effect and denominator effect can be further decomposed by decomposing 
= (2) 1) a2) «(1) 
the movement estimates for the numerator (A. -A, ) and the denominator (B - —B, ) into 


the various decomposition components (3.5 - 3.17). 

4. Decomposition "Replicate Variance" Methodology 

4.1 The Bootstrap Variance Estimator 

The bootstrap procedure for a stratified random sample selected without replacement (Shao 
and Tu 1995) is to: (i) Select simple random samples s,, without replacement (SRSWOR) of 


mn = [ny/2] units from the original sample of 1; units, independently within each stratum h. 
Let 6,4; equal 1 if ies,, and zero otherwise. Calculate the bootstrap sampling weights: 


=W; (1 = fi) (1-fi)mn n . 
=w(1 cm + [Sate (4.1) 
h 


where fj, = N is the sampling fraction within stratum h. Calculate the bootstrap estimator of 
Yo 


ve =2 WriQriYi (4.2) 


(ii) Independently replicate step (i) a large number of times, R, and calculate the bootstrap 
estimates, Vic, You ay Ves 


(iii) The bootstrap variance estimator is given by the Monte Carlo approximation: 


a iG 
x & oe ?.] 


Ro] (4.3) 


4.2 The Bootstrap Variance Estimator for Movement Estimates 


The bootstrap procedure for estimating the variance of movement estimates for stratified 
random samples selected without replacement is to: (i) Divide the units selected in the 


survey into the three categories: s\) = s\? qs}, st =s\? as) and s?° =s\? ns. (ii) Within 
each of the three categories select simple random samples 3), gil?) and 52) without 
replacement (SRSWOR) of mo =[n/2], mY = [n?/2] and m2” = [n?/2] units from the 


original samples of n©, ni) and ne?) units, independently within each stratum h. Let 69 : 


file) and joe equal 1 if ies, ies 4” and ies 2”, and zero otherwise. Calculate the bootstrap 


sampling weights: 


Cs ) Ym) Gee ) yy) n® 
of =" {1- fi), + fedmy? my 


(c s¢ 7a (c) 
TAR TSE OLS 72: a Oks Oni , if 1es 
(nj? —m’) (ni — my’) mo | ' 
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x oof | Carel —fP)- nO {oak 


= n(n) _ m) 
n\?( (1 ap f(y fe )) (tc) (te) 
(t) hh h JM Mh (te) (te) 
a [| nO(n® —m®) nor rifies,” (4-4) 
h Wy, — Mm), Mh 
(1), (2) 
ny nyoN? ; 
where f,? = 44 — _. Calculate the bootstrap movement estimator of Y"” 
h (c) 7D pg * 
Nh Nn INA 
vm) (2) 1), 
Yoo = ee Ww, ou yi 2, x wi ge Yi (4.5) 


(2 


(iii) Independently replicate step (ii) a large number of times, R, and calculate the bootstrap 
A(m 


n\n 
estimates, Y7. , Yor po VRe 


(iv) The bootstrap variance estimator is given by the Monte Carlo approximation: 


— 2 
alm alm 
5 [9-7 ) 


Var(Y.) = Roi (4.6) 
a(m) 
ain) & Yee 
where Y, = R 


4.3 The Bootstrap Variance Estimator for Decomposition Components 


The bootstrap movement estimator (4.5) can be written as the difference between a product 
of four factors: 


<(m) (2)~(2) = NMG (1)—0) — 
Lie =) Ni PincS aa them Ni P rnc rnd a CD) 


where po”, ae and y 7 : represent the bootstrap estimates of p p. i. one and a”: 


a(t) _ esi 
rhe (t) 


w” (t), (t) 
2. Wy; Bri Yi 


—(t) ies Shc 
8 rhe y w? y? 


(t) 
i he 


s(t) _ ies\!) 
rhe (t) 
x WwW, 
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Furthermore, the decomposition components (3.5 - 3.17) can be written as a function of Nv, 


a(t) <(t) —(t), 
Phes She and Vue: 


op) ey Ral) oe BOY AO 20 aay) 
Y, = y' (Ni; Ni o Pict Phe Sher her Vics Une) (4.8) 


Therefore, the bootstrap procedure for estimating the variance of decomposition 
i mee) 
components is to: (i) Calculate the bootstrap estimator of Y1., Yo. ,..., Yre, where: 


a (*) 
a7) 72) s@) 22) <M =@) —(1) —(@2) 
Yrc =W (Ni Ni VP thc E het Sricr B ret Uther Y ric) (4.9) 


(ii) The bootstrap variance estimator of the decomposition components is given by the 
Monte Carlo approximation: 


(4.10) 


QUESTION 5: Is it appropriate to use the bootstrap replicate weights to estimate the 
variance of the decomposition components? 


5. Simulation Study 


A Monte Carlo simulation study was undertaken to measure the bias and variability of the 
decomposition components and the bootstrap variance estimator of the decomposition 
components. In the Monte carlo study, an artificial population U of 12,000 units was split 


into three categories: U = U™ 4 U of 9,000 units, UW" = UY 9 U® of 1,000 units and 


u? =U > U® of 2,000 units. The characteristics of the population was generated using 
the following models: 


yP =11x|xi+y?| 
where x; = abs(C;), €; = N(O, 22.5) if ie, 6; = N(0, 25) if ic, C; = N(0, 27.5) if ic U@, 
yw? =N(0, /x;) and yw? = N(0, /x7). 


All units in the population were randomly assigned to one of two stratum groups, and then 
within these two stratum groups units were assigned to one of five strata, using stratum 
boundaries 0-9, 10-19, 20-29, 30-49, 50+ based on Zz) 7 int|x; +é& | at the first time period and 
z?) = int|x; + €| at the second time period, where €{" = N(0, /x; ) and EP) = N(O, /%i ). 

All units in the population were randomly assigned to one of two domains of interest. The 


simulation population characteristics of the variable of interest are presented in Table 3. 
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Table 3: Simulation Population Characteristics of the Variable of Interest 


Time 1 Time 2 Movement 
Time 1 Time 2 Movement Common Common Common 
Domain Population Population Population | Population’ | Population’ | Population’ 
1 97,847 122,950 25,103 55,344 61,099 5,756 
2 100,478 125,135 24,657 56,879 62,745 5,866 
Time 1 Time 2 Movement Time 1 Time 2 Movement 
Death Birth Birth Death Flip In Flip Out Stratum Flip 
Domain Population Population Population Population Population Population 
1 8,487 24,484 15,997 34,106 37,367 3,351 
2 9,263 24,689 15,426 34,336 37,701 3,365 


*The common population is defined as those units in the same stratum at both time periods. 


A total of 10,000 independent simulated stratified samples were selected from the 


population. In the first time period, stratified simple random samples without replacement 


of 200 units, with equal allocation of nv 


time period, U". In the second time period, a stratified synchronised random samples of 


220 units, with equal allocation of ni?) = 22 and no = 14, were selected from the population 


at the second time period, U. 


= 20, were selected from the population at the first 


The estimator used this study was the generalised regression estimator (2.2), with X equal to 


the population totals for the two auxiliary variables x, = {x1;,x2;} at the stratum group level, 


and cio? =], 


The relative biases (RB) and the relative root mean squared errors (RRMSE) of the 
decomposition components were estimated by: 


RB(Y.”) = si a9 000 a4 ce y)| 


RRMSE(Y,”) = a / 70-000 os (V2 -re) 


where Y, are the decomposition components (3.5 - 3.17) under the k" simulation sample 
and Y‘") are the decomposition components (3.5 - 3.17) under a complete census of the 
population at the two time periods. The estimated RB and RRMSE of the decomposition 
components for the variable of interest are presented in Table 4. 
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Table 4: Estimated Relative Biases and Relative Root Mean Squared Errors 
of the Decomposition Components 


Common Net Net Weight 
Sample Birth Stratum Net Domain Adjust 
Unit Death Flip Rotations Change Change 
Movement Effect Effect Effect Effect Effect Effect 
Domain Y. gone yo yo ahh aka ec 
Population Values 
1 25,103 10,870 13,531 235 0 468 0 
2 24,657 10,894 14,260 -11 0 -485 0 
Average Simulation Estimates 
1 24,897 10,777 13,794 -492 19 338 461 
2 25,011 10,884 14,568 -621 2 -325 505 
Estimated Relative Biases (%) 
1 -0.8 -0.4 1.0 -2.9 0.1 -0.5 1.8 
2 14 0.0 1.2 -2.5 0.0 0.6 2.0 
Estimated Relative Root Mean Squared Errors (%) 
1 29.3 6.0 9.4 8.6 5.9 26.6 7.5 
2 29.9 6.1 9.7 8.6 6.0 27.0 8.1 
The results of the simulation study found that the estimated relative biases of the 
decomposition components are quite small. These results indicate that the expected value of 


the common sample unit effect (3.11) is approximately equal to the common population unit 


effect, but it is not equal to the difference in the population values of the common 
population units: 


a (muc) 


E(Y )e yimuc) 4 yc) __ ye) 


Furthermore, the expected value of the net birth-death effect (3.10 plus 3.13 plus 3.16) and 
the net stratum flip effect (3.9 plus 3.14 plus 3.17) are approximately equal to the true 
population net birth-death effect and net stratum flip effect, but they are not equal to the 
difference in the population values of their respective population units: 


a(mfg) — ~(mub) (mud) 


EY +¥ + Yo) & Vina) 4 youd) 4 yonud) 4 y2e) _ ye) 


EY 4 sph ts ake mw Yin 4. ylonfi) 4 ynufo) 4 yf) — y fo) 


QUESTION 6: Is it necessary for the expected value of the common sample unit effect, 
the net birth-death effect, and the net stratum flip effect to be equal to the differences in 
the population values of their respective population units? 


The RB and RRMSE for the bootstrap standard error estimator of the decomposition 
components were estimated by: 
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1 
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y (se(¥)-seW)) | 


k=1 
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—_|-1_y 
SEQ) | 10, 000 


10,000 
1 


k=1 


(se(%,.)- SEQ )) 


a(x) “i 
where SE te ) are the bootstrap standard errors of the decomposition components under 


the k” simulation sample and SE(Y™) = | 


1 


10,000 7 (, 
10, 000 & ( 


) 2 
- ye) is the estimated true 


standard error of the decomposition components. The estimated RB and RRMSE for the 
bootstrap standard error estimator of the decomposition components for the variable of 
interest are presented in Table 5. 


Table 5: Estimated Relative Biases and Relative Root Mean Squared Errors for the 
Bootstrap Standard Error of the Decomposition Components 


Common Net Net Weight 
Sample Birth Stratum Net Domain Adjust 
Unit Death Flip Rotations Change Change 
Movement Effect Effect Effect Effect Effect Effect 
Donain Y. . be o nen gin a (nr) s fn) S bre) 
Estimated Population Values 
1 7,351 1,507 2,337 2,023 1,489 6,665 1,827 
2 7,358 1,516 2,375 2,023 1,483 6,649 1,933 
Average Simulation Estimates 
1 7,488 1,516 2,466 2,249 1,632 6,797 2,259 
2 7,498 1,531 2,494 2,225 1,620 6,784 2,357 
Estimated Relative Biases (%) 
1 1.9 0.6 5.6 11.2 9.6 2.0 23.6 
2 1.9 1.0 5.0 10.0 9.2 2.0 21.9 
Estimated Relative Root Mean Squared Errors (%) 
1 8.8 11.4 16.6 23.6 23.2 8.6 30.0 
2 8.8 12.1 16.9 24.1 22.9 8.7 28.2 


The results of the simulation study found that the estimated relative biases for the bootstrap 
variance estimator of the decomposition components are reasonably small, with the 
exception of the weight change adjustment effect. 


QUESTION 7: Is the size of the estimated relative biases for the bootstrap variance 
estimator of the weight change adjustment effect problematic? 
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6. Using the Decomposition Methodology to Identify Irregularities 


The decomposition methodology has dual purposes. Firstly, it can be used as a tool to better 
understand the key factors that are driving the movements in the estimates. It has the 
potential to identify problems which might influence the movements in the estimates, such 
as unexpected levels of birthing or deathing of units in the population of interest, unusual 
combinations of units rotating into and out of the sample, and errors in the auxiliary 
benchmark values or auxiliary benchmark totals. 


Secondly, it can be used as a macro editing tool to assist in the output editing of sample 
surveys. It has the potential to identify suspicious movements in the aggregate level 
estimates (using the decomposition components presented in Section 3.1 and 3.2), and then 
drill down to the micro level data to identify units with the largest impact on the movement 
estimates (using the unit level contributions to the decomposition components presented in 
Section 6.1). The suspicious movements in the aggregate level estimates can be investigated 
by looking at the size of decomposition components relative to the standard errors on the 
decomposition components as well as looking at the size of the decomposition components 
compared to the time series of the decomposition components. 


6.1 Unit Contributions to Decomposition Components 


The unit contributions to the movement estimates will have contributions from each of the 
four factors. The contribution of unit i to the frame change effect (3.5) is give by: 


ln) (2) ~@) ~@) —() —(@) ¢ +7 7(fi) 7 7b) 
ye = 9 PL Pe Be Bente) _ if icUf”, Uj 
gin _ a) g@) 50) 52) —4) <Q) Bees (aay oa 218) 
Yor =—OPrer Pe 8ne+ She Vacs Vie ) _ if icf”, Ui 


The contribution to unit i to the domain change effect (3.6) is given by: 


a2) 2) -@)_ -(1) 


oom) | Phe ~Phe_ , Phe ~Phe () py2) g@ 52) 5M 50) (0) 
Ya ond 2n?) Q(Ni, Nr Sno She Ynes Vic ) if t 1ESp, : 
im) | Bye — eal @) =@) =) =() 
Mn Phe Fhe_ (1) —(1) =(2) —(1) -@ sce (ri) off) o(b) 
ne rk 7) Q(Ni, /Ny ner 8ner Vic Vie ) , if iesy”, Si", 8} 
2nNje 
om) [Bre ~ Pye | @) =@) =@) =@) 
a(mu he Pn (1) ah 20) (1) 0 (ro) 3?) (d) 
Ya . oe : Q(Ni, /Ny She ’ She Ve Vic ) ,ift 1S), ‘ : Sh 
he J 


The contribution of unit i to the weight adjustment change effect (3.7) is given by: 


r—(2) —(1) —(2) —(1) 
ae Sine ~8ne , She ~8he | (N,N N® g® G2 ZO 2@ 


aes (c) 
h Phos Phe Vine Vne ) , if ies; 


2 
2n\ 2n? 
in [32 hig 
gnu) | oh hi (1) AV ) sf) -Q@) —G) =) (ri) of) .@) 
co ae l (Nis Ni, Paci bet Vai Mae) , if 1es;;" Sh Sh 
ne 
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=(22) =) 
ol) | Sine ~ She 1 ) a) a) =) wee ) (a) 
Ya 7 Oe (Ny, i Ny Pde ve ) : if ies, sf”, s} 
he 


The contribution of unit i to the average value change effect (3.8) is given by: 


cm) | Ye Yi | aaa) agl2) a) 2 <0) <0) ee 
ci () 7) ~(Ni Ni Phot Phe She Suc) if 1€8;, 
Nhe Nhe 
(2) —(2) 7 
clu) | Wri Yn () y@ g@ ge) 50) «cs o(ri) aff) o(b) 
Yi = 0) oz @) p(Ni,,N reNh tPhe+ Phe 1 She re ) U if i€s},, Sif Sh 
Nhe Nhe J 
ae tye al ett 
yams | She Shi (1) x72) =@) »@) =<G) =@) +c _(ro) o(fo) o(d) 
ci Ho: yO ON Ni Phe Phe Shes She ) , ifies;” sexs 
he he 4 


6.2 Ability of the Decomposition Methodology to Identify Irregularities 
(i) Example 1: Unusual Rotation Combinations 


A single simulation sample from the simulation study was used to examine the impact of 
unusual combinations of units rotating into and out of the sample. Various combinations of 
the smallest and largest units were rotated into and out of the sample at the second time 
period. The estimates and standard errors on the decomposition components of using these 
unusual rotation combinations are presented in Table 6. 


Table 6: Impact on Decomposition Components of Unusual Rotation Combinations 


Net Net Weight 
Common Birth Stratum Net Domain | Adjust 
Time 1 Time 2 | Movement | Sample Death Flip Rotations | Change | Change 
Rotation Estimate | Estimate | Estimate Effect Effect Effect Effect Effect Effect 
Combinations ran 7?) al gine) ge qn gh) ge i 
No 103,302} 129,965 26,664*| 12,004") 14,533* -206 474 -468 307 
adjustments (7,197) (1,467) (2,580) (2,679) (1,524) (6,539) (2,213) 
Smallest units} 103,302) 133,915 30,614") = 12,745") = 13,598* -133 4,383* 2,358 -2,337 
rotated out (7,176) (1,573) (2,531) (2,529) (1,847) (6,991) (2,167) 
Largest units 103,302] 125,404 22,102") 11,596") = 14,770* -470 -6,616* -246 3,068 
rotated out (8,227) (1,526) (2,599) (2,368) (2,323) (6,446) (2,302) 
Smallest units} 103,302) 122,518 19,216*} 11,955*| =14,109* -81 -7,259* -2,739 3,230 
rotated in (6,965) (1,473) (2,636) (2,685) (2,902) (6,295) (2,189) 
Largest units 103,302} 125,503 22,202") 11,407*) 14,089* -104 7,598* -5,787 -5,002 
rotated in (8,188) (1,424) (2,490) (2,594) (3,287) (6,748) (2,648) 


* Significantly different from zero at the 95% confidence level. 
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The unusual rotation combinations had substantial impacts on the net rotations effect. Since 
the net rotations effect were significantly different from zero at the 95% confidence level, 
further investigations would be undertaken in order to identify the underlying reasons 
behind the large net rotations effects. The unusual rotation combinations also had an large 
impact on the weight adjustment change effect. The larger weight adjustment change effects 
occurred because units with the very small and very large auxiliary benchmark values were 
rotated into and out of the sample at the second time period, since the reported values are 
highly correlated with benchmark variables. 


(ii) Example 2: Invalid Benchmark Totals 


A single simulation sample from the simulation study was used to examine the impact of 
errors in the auxiliary benchmark totals. The auxiliary benchmark totals at the second time 
period were modified by constant factors ranging between 0.95 to 1.05. The estimates and 
standard errors on the decomposition components of using these invalid auxiliary 
benchmark totals are presented in Table 7. 


Table 7: Impact on Decomposition Components of Invalid Benchmark Totals 


Net Net Weight 
Benchmark Common Birth Stratum Net Domain Adjust 
Total Time 1 Time 2 | Movement) Sample Death Flip Rotations | Change Change 
Adjust Estimate | Estimate | Estimate Effect Effect Effect Effect Effect Effect 
Factor y” oe 7” ye i ve i ce 7 
1.00 100,141 125,212 25,072* 11,805* 15,780* -120 -317 -1,811 -267 
(8,400) (1,729) (2,516) (2,271) (1,555) (7,696) (2,088) 
0.95 100,141 92,474 -7,666 10,269* 13,598* -120 -258 -1,722|  -29,434* 
(6,910) (1,541) (2,198) (2,047) (1,349) (6,619) (2,839) 
0.99 100,141} 117,980 17,840* 11,466* 15,298* -120 -304 -1,791 -6,710* 
(8,000) (1,687) (2,445) (2,222) (1,509) (7,457) (2,049) 
1.01 100,141} 132,815 32,674* 12,162* 16,287* -120 -331 -1,831 6,507* 
(8,853) (1,774) (2,590) (2,325) (1,603) (7,947) (2,266) 
1.05 100,141} 167,242 67,102* 13,778* 18,582* -120 -393 -1,924 37,179* 
(11,214) (1,978) (2,932) (2,567) (1,824) (9,089) (4,039) 


* Significantly different from zero at the 95% confidence level. 


Errors in the auxiliary benchmark totals as little as one percent had substantial impacts on 
the weight adjustment change effect. Since the weight adjustment change effect were 
significantly different from zero at the 95% confidence level, further investigations would be 
undertaken in order to identify the underlying reasons behind the large weight adjustment 
change effect. 
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(ii) Example 3: Extreme Reported Values 


A single simulation sample from the simulation study was used to examine the impact of 
extreme reported values. A total of four units; one common sample unit (#01011), one 
sample rotation in (#02609), one sample flip in (#01500) and one sample birth (#01080); were 
modified by adding 100 to their reported values at the second time period. The estimates 
and standard errors on the decomposition components of using these four units with 
extreme reported values are presented in Table 8. 


Table 8: Impact on Decomposition Components of Extreme Reported Values 


Net Net Weight 
Common Birth Stratum Net Domain Adjust 
Time 1 Time 2 | Movement | Sample Death Flip Rotations | Change | Change 
Extreme Estimate | Estimate | Estimate Effect Effect Effect Effect Effect Effect 
Values g (2) om) oAmuc) (mba) gin rai (md) ye 
No extreme 96,651) 123,561 26,910* 10,570* 17,905* -906 -475 280 -463 
values (6,939) (1,516) (2,675) (2,618) (1,114) (6,163) (2,318) 
Extreme 96,651) 146,833 50,182*} = 20,379* 17,856* 4,482 4,808 1,925 733 
values (14,536) (9,976) (4,045) (6,914) (4,850) (6,553) (2,593) 


* Significantly different from zero at the 95% confidence level. 


Although these four units with extreme reported values did not result in the net stratum flip 
effect or the net rotations effect to be significantly different from zero at the 95% confidence 
level, further investigations would be undertaken in order to identify the underlying reasons 
behind the unexpectedly large movement estimate. The ten largest unit contributions to the 
movement estimates using these four extreme reported values are presented in Tables 9 and 


10. 
Table 9: Largest Unit Contributions to Movement Estimates 
Without Extreme Values 
Weight 
Frame Domain | Adjustment Unit 
Time 1 Time 2 | Decomposition Change Change Change Value Total 

Unit Value Value Effect Contribution | Contribution | Contribution | Contribution | Contribution 
A N.a. 13.6/Flip In 18.3 115.3 -11.8 -1,392.9 -1,271.1 
B 28.2 n.a.|Flip Out -5.4 -38.0 14.3 -1,124.9 -1,154.0 
Cc N.a. 64.0|Birth 18.3 115.3 -11.8 765.8 887.6 
D 27.3 n.a.|Flip Out -18.3 144.1 -14.7 774.2 885.2 
E 53.3 60.3}Common Unit 0.0 259.4 -26.5 596.9 829.8 
F N.a. 24.4|Birth 5.4 -38.0 14.3 803.6 785.2 
G 59.2 64.4|Common Unit 0.0 259.4 -26.5 448.4 681.2 
H 46.6 n.a.|Rotations Out -23.1 -89.7 -23.9 -533.9 -670.6 
I 38.8 43.8|Common 0.0 259.4 -26.5 430.3 663.2 

J 56.2 n.a.|Flip Out -18.3 144.1 -14.7 -769.4 -658.3 
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Table 10: Largest Unit Contributions to Movement Estimates 


With Extreme Values 
Weight 
Frame Domain | Adjustment Unit 
Time 1 Time 2 | Decomposition Change Change Change Value Total 
Unit Value Value Effect Contribution | Contribution | Contribution | Contribution | Contribution 
K 5.5 105.5|Common Unit 0.0 274.1 181.1 9,810.9 10,266.1 
L n.a. 125.5]Flip In 10.2 6.9 9.6 6,566.8 6,593.5 
M n.a. 100.6|Rotations In 6.5 112.1 74.1 4,968.6 5,161.3 
N n.a. 166.2|Birth 32.9 -45.4 -37.3 1,510.9 1,461.2 
A n.a. 13.5/Flip In 18.3 115.3 -11.8 -1,392.9 -1,271.1 
B 28.2 n.a.|Flip Out -5.4 -38.0 143 -1,124.9 -1,154.0 
Cc n.a. 64.0|Birth 18.3 115.3 -11.8 765.8 887.6 
D 27.3 n.a.|Flip Out -18.3 144.1 -14.7 774.2 885.2 
E 53.3 60.3}Common Unit 0.0 259.4 -26.5 596.9 829.8 
O n.a. 1.6|Birth 6.5 112.1 74.1 -1,003.3 -810.5 


The four units with extreme reported values were the four largest contributors to the 
movement estimates, and since the contribution of three of these units were much greater 
than the other units, these units would be given the highest priority in the output editing of 
the survey. 


7. Conclusion 


The proposed decomposition of the movement estimates into primary and secondary factors 
is just one of many possible approaches to the decomposition of the movement estimates. 
This decomposition methodology is consistent with the generalised regression estimator 
currently used in ABS business surveys, and the existing bootstrap variance estimation 
methodology can be used to produce approximate variances for the various decomposition 
components. 


This decomposition methodology can be used as a tool to better understand the key factors 
that are driving the movements in the estimates and as a macro editing tool to assist in the 
output editing of sample surveys. It has the ability to identify suspicious movements in the 
aggregate level estimates and then to drill down to the micro level data to identify units 
with the largest impact on the movement estimates. 


This decomposition methodology has been implemented in a number of ABS Business 
Surveys, although the variances for the decomposition components have yet to be 
implemented. Although the decomposition methodology has proven useful to users, the 
users have had some difficulties with the interpretation of the decomposition components 
and understanding the significance of the size of the decomposition components. It is 
envisaged that the availability of variances for the decomposition components, longer time 
series of the decomposition components, and training on the interpretation of the 
decomposition components might alleviate this limitation. Another limitation of the 
decomposition methodology is the lack of a link to the final published seasonally adjusted 
movement estimates. 
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Appendix: Descriptions of the Movement Decomposition Components 


Decomposition 
Component 
Frame Change 
Effect (3.5) 


Description 
Measures the impact on 
the movement estimates of 
changes in the number of 
units on the survey frame 
between the two time 
periods. 


Positive Effect 


Negative Effect 


Frame Growth 
Effect (3.9) 


Measures the impact on 
the movement estimates of 
the net difference between 
the number of units 
birthed and deathed on 
the survey frame between 
the two time periods. 


Larger number of units 
birthed on the survey frame 
than deathed on the survey 
frame, or units deathed on 


the survey frame occur within 


strata with a lower average 
value of units, or units 
birthed on the survey frame 
occur within strata with a 


higher average value of units. 


Smaller number of units 
birthed on the survey frame 
than deathed on the survey 
frame, or units deathed on 
the survey frame occur within 
strata with a higher average 
value of units, or units 
birthed on the survey frame 
occur within strata with a 
lower average value of units. 


Frame Stratum Flip 
Effect (3.10) 


Measures the impact on 
the movement estimates of 
the net difference between 
the number of units on the 
survey frame that flipped 
into and out of strata 
between the two time 
periods. 


Units have flipped out of 
strata with a lower average 
value of units and flipped 
into strata with a higher 
average value of units. 


Units have flipped out of 
strata with a higher average 
value of units and flipped 
into strata with a lower 
average value of units. 


Domain Change 


Measures the impact on 


Higher estimated proportion 


Lower estimated proportion 


Effect (3.6) the movement estimates of | of units within domain of of units within domain of 
changes in the estimated _| interest at the second time interest at the second time 
proportion of units within | period than at the first time period than at the first time 
the domain of interest period. period. 
between the two time 
periods. 

Weight Adjustment | Measures the impact on Higher weighted average Lower weighted average 

Change Effect (3.7) | the movement estimates of | g-weights at the second time | g-weights at the second time 


changes in the weighted 
average g-weight between 
the two time periods. 


period than at the first time 
period. 


period than at the first time 
period. 


Average Value 
Change Effect (3.8) 


Measures the impact on 
the movement estimates of 
changes in the average 
value of units of the 
variable of interest 
between the two time 
periods. 


Common Sample 
Unit Effect (3.11) 


Measures the impact on 
the movement estimates of 
differences in the average 
value of common sample 
units between the two 
time periods. 


Higher average value of 
common sample units at the 
second time period than at 
the first time period 


Lower average value of 
common sample units at the 
second time period than at 
the first time period 
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Decomposition 
Component 


Description 


Positive Effect 


Negative Effect 


Sample Rotations 


Measures the impact on 


Higher average value of 


Low average value of sample 


In Effect (3.12) the movement estimates of | sample rotations in than rotations in than common 
differences in the average |common sample units at the /sample units at the second 
value of sample rotations |second time period. time period. 
in and common sample 
units at the second time 
period. 

Sample Stratum Flip | Measures the impact on Higher average value of Low average value of sample 

In Effect (3.13) the movement estimates of | sample flips in than common | flips in than common sample 
differences in the average |sample units at the second units at the second time 
value of sample that time period. period. 
flipped into strata and 
common sample units at 
the second time period. 

Sample Birth Measures the impact on Higher average value of Low average value of sample 

Effect (3.14) the movement estimates of | sample births than common ___|births than common sample 


differences in the average 
value of sample births and 
common sample units at 
the second time period. 


sample units at the second 
time period. 


units at the second time 
period. 


Sample Rotations 


Measures the impact on 


Lower average value of 


Higher average value of 


Out Effect (3.15) the movement estimates of | sample rotations out than sample rotations out than 
differences in the average |common sample units at the |common sample units at the 
value of sample rotations | first time period. first time period. 
out and common sample 
units at the first time 
period. 

Sample Stratum Flip | Measures the impact on Lower average value of Higher average value of 

Out Effect (3.16) the movement estimates of | sample flips out than sample flips out than 
differences in the average |common sample units at the |common sample units at the 
value of sample that first time period. first time period. 
flipped out of strata and 
common sample units at 
the first time period. 

Sample Death Measures the impact on Lower average value of Higher average value of 

Effect (3.17) the movement estimates of |sample deaths than common _|sample deaths than common 


differences in the average 
value of sample deaths 
and common sample units 


at the first time period. 


sample units at the first time 
period. 


sample units at the first time 
period. 
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